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ABSTRACT. We employ a modified cascade hydrodynamics code to simulate 
the phase transition of an expanding quark-gluon plasma and the passage of 
a charmed particle through it. When inside the plasma droplets, the charmed 
quark experiences drag and diffusion forces. When outside the plasma, the 
quark travels as a D meson and experiences collisions with pions. Additional 
energy transfer takes place when the quark enters or leaves a droplet. We find 
that the transverse momentum of D mesons provides a rough thermometer 
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1 Introduction 



Charmed particles are promising candidates for probing the early stages of 
ultrarelativistic heavy-ion collisions. Because of their large mass, they can 
be created only in the initial hard collisions between the nuclear constituents 
[p]]. Most charmed quarks created in hadronic collisions emerge with non- 
relativistic velocities 0], and hence in nuclear collisions they will remain in 
the collision region for a long time. Finally, the passage of charmed quarks 
through the plasma is a problem amenable to theoretical analysis [|, |5], 
Thus the observation of charmed particles, and in particular the measure- 
ment of their transverse momentum distribution, should yield information 
about the hot environment present early in the system's evolution. 

Consider a high-energy nuclear collision which leads to creation in the cen- 
tral rapidity region of both a charmed quark pair and a quark-gluon plasma. 
We adopt Bjorken's scaling hydrodynamics ansatz to describe the plasma 
[0, The initial creation of the cc pair takes place in a time on the or- 
der of l/m c ~ 0.1 fm/c. Within a few tenths of fm/c the quarks will find 
themselves in an expanding quark-gluon plasma near local equilibrium 
We begin following the progress of each quark at this time, r = r , when 
the temperature is T . At this point, we assume, each quark begins a pro- 
cess of scattering and diffusion in the plasma which causes a loss of initial 
momentum and relaxation towards the (changing) thermal velocity. 

The plasma expands and cools according to T = To • (r/ro) -1 ^ 3 , until at 
r = t p the phase transition begins, the temperature having reached T = T*, 
the transition temperature. The "mixed phase" forms as follows. If there 
is a genuine first-order transition, the plasma supercools slightly, and then 
bubbles of the low-temperature phase are nucleated by fluctuations or by 
inhomogeneities. These bubbles grow, meet, and coalesce, and eventually 
we have a hadron gas containing mostly pions and also droplets of plasma 
which are evaporating. The phase transition is complete when the last of the 
droplets evaporate into hadrons, at r = th- The enormous entropy of the 
plasma which must be converted to hadrons dictates that is more than 
ten times r p : The mixed phase is the longest stage in the evolution of the 
system. 

At the beginning of the mixed phase each charmed quark is in plasma, 
but it may emerge into a region of hadron gas, and then it might be reab- 



sorbed into plasma only to emerge again [10|. Emergence into the hadron 
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gas requires hadronization, the binding of the charmed quark with a light 
antiquark to form a meson; reabsorption involves stripping off the antiquark 
by collisions with the plasma. We will discuss the processes of hadronization 
and ionization in detail below; suffice it for now to observe that hadroniza- 
tion involves a loss of kinetic energy by the quark, enough to supply at least 
the difference mo — m c — 350 MeV. If the quark does not have this kinetic 
energy to lose, it may be trapped in the plasma. It could emerge later if it 
gains enough energy through diffusion; otherwise it is trapped in the plasma 
until the plasma finally evaporates. Of course, it must hadronize at least 
once. Whenever the quark/meson is in a region of hadron gas, it is liable 
to scatter off the hadrons, mostly pions, and thus to undergo yet more drag 
and diffusion. 

The situation is not too different if there is no real phase transition. The 
high- and low-temperature regimes are still distinguished by very different 
entropy densities. The "mixed phase" in this case would be a mixture of 
regions with temperatures slightly above and slightly below the crossover 
temperature at which the entropy changes steeply. The notion of bound- 
aries between droplets of the pure phases should still be valid, as should the 
estimates of r p and th- 

For simulation of the mixed phase we adopt the cascade hydrodynamics 
model [T]J of Bertsch et al. In this model, when the temperature reaches 
T* the plasma simply breaks up into droplets. As the system continues 
its longitudinal expansion the droplets become ever more rarefied, all the 
while radiating and reabsorbing pions. Finally the droplets all evaporate, 
completing the phase transition. We have supplemented the cascade code 
with the routines needed to follow the progress of a charmed quark/meson 
through the system. 

In the next section we describe the various assumptions of our model 
and how they are implemented in the simulation program [12]. Proceeding 
step-by-step through the collision process, we discuss the initial conditions 
given the plasma and the charmed quark diffusing therein; the Langevin 
process by which this diffusion takes place; and the cascade hydrodynamics 
by which the plasma droplets and pions are simulated in the mixed phase^. 
We then present the model which is central to the physics presented here, 
that of charm hadronization and re-ionization [10], and finally discuss the 

3 This is essentially unchanged from the algorithm originally used by Bertsch et al. |0| 
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Dtt scattering which occurs when the charmed quark is outside the plasma. 

In Section 3 we present our results, which are mainly predictions for the 
final p± of a charmed meson born of a c quark created with initial p± = p±Q. 
We plot this quantity for a range of assumptions concerning the free param- 
eters of the simulation. Our main qualitative result may be simply stated: 
D mesons do provide a rough thermometer of the phase transition. This is 
because a large percentage of the charmed quarks are trapped in the plasma 
until it breaks up; moreover, those which escape the plasma droplets early 
will experience many collisions in the pion gas which has not yet rarefied. 
Unfortunately, the temperature scale of the thermometer contains uncer- 
tainties due to the dynamical assumptions we make, chiefly the initial size of 
the plasma droplets. The measurement of temperature is also degraded by 
transverse flow of the droplets and pions. 



2 Ingredients of the simulation 
2.1 Initial conditions 

The initial conditions we specify are of two classes. One describes the cylinder 
of plasma at proper time r , while the other gives the initial momentum and 
location of each charmed quark. 

The plasma is assumed to occupy a cylinder of fixed cross section R, equal 
to the radius of the colliding nuclei. We consider only central UU collisions, 
and hence R ~ 7 fm. The initial time tq is the time at which the plasma is 
taken to be fully formed. We can estimate it from the flux tube model P~3" 
according to [Q r ~ (1 fm) • A -1 / 6 ~ 0.5 fm. 



We fix the initial temperature in the usual way [|T5[ by setting the eventual 
pion multiplicity dN n /dy. This is related to the entropy density at freezeout 
by dS/dy = 3.6 dN n /dy. The cascade simulation predicts an entropy increase 



of 20-30% during the phase transition JTT]]; if we take this into account, 
while neglecting the entropy created through dissipation at other stages of 
the expansion, then the entropy density dS/dy at To is (1.2) -1 times that at 
freezeout. Setting 

dy \ J 8J 45 w 
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in the quark-gluon plasma, we determine the temperature To at To. (For 
simplicity we take Nf = 2.) We present results for an initial temperature of 
300 MeV, corresponding to dN/dy = 1400. 

The hard collision of the incident nuclei takes place at t = z — 0, and 
hence the charmed quarks are created on trajectories emanating from this 
surface and thus satisfying z = v$t for some v%. We follow one quark at a 
time, and so the longitudinal-boost invariance of the fluid allows us to work 
in the frame where vg = 0. By choosing the x axis along the quark's initial 
position we end up, with no loss of generality, withF] x = (xo, 0, 0) and 
Po = (pS,Po,0). 

The radius Xq at which the heavy quark is created is taken from a distri- 
bution proportional to the number of nucleon-nucleon collisions taking place 
at that radius, 

p{xq) dx oc (^R 2 — x^j 2tcx dx . (2) 

Azimuthal symmetry of the parton-parton collision gives the direction of 
Po in the xy plane a uniform distribution. The magnitude p is an input 
parameter. 

Charmed quarks are of course created in pairs. We neglect entirely the 
interaction between the quark and antiquark, and any influence they may 
have on each other as they diffuse through the plasma and mixed phases. 
The reason is that the two are created with opposite transverse momenta 
(and different rapidities) which carry them apart quickly. Since r p is (under 
our assumptions) at least 1.7 fm/c, the quarks can be 3 fm apart by the onset 
of the phase transition; since the droplets in the mixed phase have initial radii 
of 1-2 fm (see below), the quark and antiquark will generally find themselves 
in different droplets.Q Their interaction will therefore be indirect at best. We 



have not studied correlations in the directions of the produced D pair [|1G 
which may offer another diagnostic probe of the phase transition. 



4 Actually for reasons of economy we put two charmed quarks into each event. We 
make sure that they have no effect on each other by giving them different initial values of 
Xq, with opposite sign, and the same po- 

5 The rare case where the cc pair are found in the same droplet might make a significant 
contribution to J/ip production. 
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2.2 Diffusion before the phase transition 

From t = m^ 1 to r = r the interaction region is in a rapidly changing 
state of particle creation and equilibration. Since this period is very short, 
it doesn't make much difference how we model it. For simplicity, we proceed 
as if it were an equilibrium plasma at T = To. We allow the charmed quark 
to diffuse in this plasma in the same way as it does after To, which we now 
describe. 

From To to t p the charmed quark is in a pure plasma. We model its 
diffusion by a non-relativistic Langevin equation, 

^ = -7(T)p + T7, (3) 

where 77 is a Gaussian noise variable, normalized such that 

(r ll {t)r 1] {t'))=a{T)5 l3 5{t-t>) . (4) 

Both the drag coefficient 7 and the momentum-space diffusion coefficient a 
depend on the local temperature. We take 7 from the result of Q , which we 
parameterize as^ 

7(T) = aT 2 , a = 2 x 1(T 6 far 1 MeV~ 2 . (5) 

(We neglect any momentum dependence in 7.) 7 and a are related by the 
fluctuation-dissipation relation, which says that in equilibrium 

2 \ « 



For temperatures where the quark in equilibrium is non-relativistic, we would 
have (p 2 ) = m c T; in the range T =150-250 MeV, however, (pf) is up to 33% 
larger than this. As a relativistic correction, we calculate a from (||) with 
(p 2 ) = 1.33 m c T. 

We apply (0), which is not Lorentz-invariant, in the rest frame of the 
plasma surrounding the quark. In the Bjorken scaling expansion we assume, 
this is a frame moving longitudinally with v z = z/t, where z and t are the 
coordinates of the quark. The local temperature is given by T = Tq(tq/t) 1 I z . 

6 We here correct an error in |ITo|l. 
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2.3 Simulation of the mixed phase 



The cascade hydrodynamics model |TTJ of the mixed phase assumes that at 
proper time t p , the beginning of the phase transition, the plasma breaks 
up into droplets. These droplets subsequently emit and absorb pions. The 
droplets' interior is always at the transition temperature T* and thus their 
energy density is fixed at 

+6iV/ DS T * 4+5 ' (7) 

where B is the bag constant. We worked with two values for the transition 
temperature, T* = 150 MeV and T* = 200 MeV, with B adjusted accord- 
ingly. Table [I] shows the values of various quantities for the two cases. 
A droplet with radius has mass equal to 

M d = ^nr 3 d e . (8) 

The initial value of r d determines how long it takes the droplets to evaporate. 
If the transition were adiabatic, the lifetime of the mixed phase would be 
determined by entropy conservation, which makes the ratio of times equal to 
the ratio of entropy densities, 

^ = ^ ~ 10 . (9) 

Tp Oh 

In the cascade model, however, the expansion is not adiabatic and the dura- 
tion of the transition, that is, the time needed for the droplets to evaporate, 
is strongly dependent on the droplets' initial radius. Thus larger initial radii 
give a longer-lived mixed phase. [] We experimented with different initial val- 
ues for r d ranging from 1 fm to 2 fm. As each droplet emits and absorbs 
pions its mass changes, and r d changes accordingly. 

The program actually steps through the time coordinate t rather than 
the proper time r. At each time t we determine the z coordinates of the two 
points on the hyperbola r = t p , where r = \J z 2 — t 2 : This is where the phase 
transition is currently occurring, and thus where droplet creation should take 

7 We find that for ro = 1.0, 1.5, and 2.0 fm the duration of the phase transition turns 
out to be 12, 22, and 32 fm/c, respectively, independent of the transition temperature. 
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Transition temperature 



150 MeV 



200 MeV 



Beginning time of transition 
Bag constant 
Plasma energy density 
Thermal pion energy 
Thermal pion momentum 
Thermal D momentum 



t p 

£ 
T 

pi 

T 

Pd 



(fm) 
(MeV) 
(GeV/fm 3 ) 
(MeV) 
(MeV) 
(MeV) 



4.0 
208 
1.0 
480 
520 
1010 



1.7 
278 
3.2 
610 
680 
1200 



Table 1: Physical quantities for the two transition temperatures considered, 
and p T D are root-mean-square momenta, p T = J (p 2 ), while E T = (E). 



place. The droplets are uniformly distributed across the cylinder, and their 
number in the £-slice at r v is fixed to conserve the entropy of the plasma at 
the breakup. A droplet's initial momentum is thermally distributed in the 
locally comoving rest frame, which has rapidity y = tanh _1 (z/£). 

From the time of its formation, each droplet interacts with the pion gas. 
By detailed balance, the pion emission rate is equal to the absorption rate 
in equilibrium |T7|. The latter is determined by assuming that a droplet 



absorbs all pions that approach within a distance d = 1 fm from its surface. 
The absorption rate in equilibrium is thus 1/4 of the pion flux times the area 
of the extended droplet surface, 



We take 



W = ^v ir -Air(r d + d) 2 . (10) 
16 

n? = 3-C(3)T 3 (11) 

and v n = 1, as appropriate for ultrarelativistic pions. Each emitted pion 
is created at a random place on the droplet's surface, with a thermally dis- 
tributed momentum (in the droplet's rest frame) directed outward. The 
emission process conserves energy and momentum exactly through the recoil 
and shrinkage of the rigid droplet. Pion absorption is simpler, consisting 
of the disappearance of any pion which ventures within a distance d of a 
droplet's surface (in its rest frame), with its four-momentum taken up by the 
droplet's recoil and growth. 

If the droplet contains less rest energy than that required to make four 
thermal pions (see Table [I]), it falls apart altogether into n thermally dis- 
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tributed pions. n might in fact be less than four, because if Md/m^ < 4 then 



n is given by the integer part of this ratio. As in |TT|, we adopt the simple 
model of first giving each pion a random, thermally-distributed momentum; 
the momenta are then shifted and rescaled to conserve the momentum and 
energy of the droplet. 

7T7T scattering is taken to be isotropic, with an isospin- and energy-depend- 
ent cross section of about 30 mb.0 Droplet-droplet scattering is neglected 
entirely, which is reasonable after the earliest times because dilution by the 
longitudinal expansion is rapid. 

2.4 Charm hadronization and re-ionization 

At the beginning of the phase transition the charmed quark is inside a droplet. 
(We restart the simulation if the random creation of droplets leaves the quark 
outside.) It continues to undergo diffusion in the plasma as it did before the 
transition time. The Langevin equation (|]) is applied in the droplet's rest 
frame via the appropriate Lorentz transformation, and every change in the 
quark's momentum is balanced by recoil of the droplet (and a change in its 
mass to conserve energy). Sooner or later, the quark strikes the wall of the 
droplet. 

Our model of hadronization is based on a picture of a quark incident on a 
static bag boundary ]18|, |19|, jSOfl. Upon striking the wall the quark stretches a 



flux tube [[B| outward which we assume to be always radially oriented. The 
tube possesses a tension (energy density) a = 0.16 GeV 2 in the droplet rest 
frame, and we take this to be the only force acting on the quark during this 
period. Thus we integrate the equation 

in the droplet frame until either the quark reenters the droplet or the flux 
tube breaks.^ 

The flux tube fissions at a rate dT/d£ per unit length. Various estimates 



for this quantity are 0.3 fm from the flux tube model |19| and anywhere 



from 0.5 to 2.6 fm 2 from comparisons of string fragmentation models with 

8 This is the weighted average at T = 200 MeV. See fll| for details. 
9 The droplet also reacts to the tension of the tube. 
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experiment We adopt two values for dT/d£ at the extremes of the range, 
0.5 and 2.5 fin" 2 . 

The tube breaks at a point V some distance £ from its base. The segment 
from V to the droplet is reabsorbed in the droplet; the remainder of the tube 
snaps back into the quark and turns it into a meson. The D emerges with 
energy (in the droplet rest frame) E^ = E c + o(L — £), where L is the total 
length of the tube. Clearly o(L — £) cannot be less than mu — E c , or there 
will not be enough energy for hadronization at all.Q 

A quark which strikes the droplet wall with energy E c < nip would be 
unable to hadronize if the droplet were static. Its reflection from the droplet 
wall, however, still involves stretching a flux tube and this process can take 
1.5 fm/c in the droplet rest frame. During this time, the droplet emits and 
absorbs pions and thus its recoil might stretch the flux tube further.^ This 
effect, which of course applies as well to situations with E c > trd, sometimes 
enables even quarks which are initially below threshold to hadronize. Again, 
the flux tube breaks at a random point, with the constraint a(L — £) > 
t(Ib — E c in the droplet 's new rest frame. 

When the tube breaks, the resulting D meson emerges immediately with 
its momentum in the same direction as that of the quark. Four-momentum 
is conserved via recoil of the droplet. In the droplet's initial rest frame, the 
kinetic energy of recoil comes from the attached segment of the flux tube and 
from a change in its rest mass, 

M' d 2 + (Ap) 2 = (M d + a£) 2 , (13) 

where Ap = p c — po is the momentum transfer from the quark/meson. 

It is possible for the D-meson to be created with momentum directed 
towards the droplet, and even for it to collide with the droplet in the next 
time step. At that point the quark is reabsorbed in an ordinary ionization 
process (see below), which means it has been effectively reflected back into 
the droplet. 

We neglect entirely an additional possibility for hadronization, namely, 
that the c quark might encounter a light antiquark as it nears the droplet 



10 We neglect the possibility of fragmentation of the quark into a D* or a Dir pair, 
because of the sharply reduced phase space for such fragmentation. 

11 A change in the droplet's radius also causes the flux tube to change its length, con- 
serving energy via the droplet's mass. 
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surface. It could thus hadronize without forming a flux tube, and possibly 
even gain energy in the process. We calculate the probability of this to 
be small, as follows. The density of light antiquarks at T = 200 MeV is 



n q = 1 fm -3 , whence the flux of outward-going antiquarks at the droplet 
surface is j q = n q /A. An antiquark must appear within a range R ~ 0.5 fm 
of the c quark during the reflection process which takes At ~ 1.5 fm/c. The 
relative probability of forming a color singlet is 1/9. Thus the probability of 
hadronization through light-quark binding is only 

P=3 q -(nR 2 )-At-l~3%. (14) 

The only remaining possibility for hadronization is that the droplet breaks 
up around the c quark, which happens when the mass of the droplet falls be- 
low the energy needed to make four thermal pions plus the energy needed for 
hadronization. In our scenario, the charmed quark picks up matter which was 
at rest in the droplet rest frame, and thus its three-momentum is conserved, 
Pd = Pc, while the energy required for creation of the D is supplied by the 
droplet. Having given this energy to the meson, the droplet then breaks up 
according to the rule given above for ordinary droplet breakup. 

If the c quark is inside a flux tube at the time the droplet breaks up, the 
flux tube's energy is merely added to that of the droplet, and the remainder 
of the breakup is handled as above. 

The inverse of hadronization, ionization, occurs when a D meson strikes 
a droplet, and is easily dealt with. We assume that the light quark is simply 
stripped off the meson, leaving a c quark with the same velocity v c = vo 
(a frame-independent statement). This means that the quark/meson loses 



energy since E c = {m c /mr))E]j. As noted in [IQ|, a D meson which undergoes 
ionization and re-hadronization in quick succession would lose energy in the 
amount 

AE = AE ion + A£ had = E D ( 1 - — ) + 

V m D J 

/ 777 \ 

> m D 1 1 ~ 370 MeV . (15) 

V m D J 

A quark with less kinetic energy than this will get trapped in the droplet. Of 
course, our simulation actually follows in detail the diffusion of the c quark 
inside the droplet so the estimate (pT5|) is only a rough guide. 
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2.5 Scattering in the hadron gas 



Whenever the charmed quark is outside a droplet it moves as a D meson. It 
can then scatter off the pions which have been emitted from the droplets. We 
use a simple model of isotropic, energy-independent, and charge-independent 
scattering, with total cross-section od-k = 9 mb. We arrive at this number via 
the additive quark model, as follows. Beginning with the (high-energy) pp 
cross section, a pp = 40 mb, we estimate that the cross section for light quark 
scattering is a qq = cr pp /9 ~ 4 mb. Similarly we begin with the ipN cross 



section f22]l , a^N = 2 mb, and estimate that the cross section for charmed 
quarks on light quarks is a cq ~ 0.3 mb <C <y qq . Thus we arrive at ajjn = 
2a qq + 2a cq ~ 9 mb. We ignore resonant scattering via the D* since the 
resonances lie at p = 39-44 MeV/c in the center-of-momentum frame, which 
is far below the typical thermal momenta (see Table [Xj) . 

(This procedure gives good results for Kir scattering. Here one begins 
with (j^tv = 13.8 mb |^3| and reaches ok-k — 13 mb, which is not far from 
the values from phase shift analysis [^4] 8 mb< ok-k < 13 mb in the energy 
range 1 GeV< ^/s < 1.5 GeV, above the K* resonance.) 



3 Numerical results 

The aim of our calculation is a prediction for the pj_ distribution of D mesons 
created in heavy ion collisions. The result naturally depends on the distri- 
bution one assumes for the initial momentum of the c quark. Rather than 
tie ourselves to specific model predictions for the latter p, we present our 
results as plots of the RMS p^ for given initial quark momentum p c ±Q in the 
range between 1 GeV and 2 GeV. 

Our main result is apparent in the last figure, which shows that p® can 
be a thermometer of the mixed phase. We begin, however, with a detailed 
analysis of the numerical results for T* = 150 MeV. In Fig. [I] we show the cal- 
culated p^ for two values of the initial droplet radius ro and of the flux-tube 
fission rate dT/d£, as discussed above. We note two features immediately: 
(1) p® varies but weakly with p c ±0 , showing the effects of drag and thermal- 
ization;0 (2) p® lies generally above its thermal value ^/§Pd = 820 MeV (see 

12 The dependence will be even weaker for p c ±0 < 1 GeV because thermalization will 
wipe out any memory of the initial momentum. 
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Figure 1: RMS transverse momentum of D meson vs. initial transverse 
momentum of c quark, for transition temperature T* = 150 MeV. The four 
curves correspond to different values of ro and dT/d£ (in fm and fm -2 , respec- 
tively): [1.0, 0.5] (dotted), [1.0, 2.5] (short dashes), [1.5, 0.5] (long dashes), 
[1.5, 2.5] (dashed-dotted). A typical statistical error bar is shown. Results 
with r = 2.0 fm (not shown) are close to, and just below, the corresponding 
points for r$ = 1.5 fm. 

Table 1), showing the effects of transverse hydrodynamic expansion. More- 
over, for any given p c ±0 , the width of the distribution of p^ (not shown) is 
consistent with a thermal distribution at T = 150 MeV. 

In comparing the four data sets in Fig. 0, we find that raising either ro 
or dT/d£ causes a drop in (for p c ±Q < 1700 MeV). We explain this by 
selecting two populations of D mesons: those which emerge from flux-tube 
fission and then escape the system — fragmentation mesons — and those which 
are trapped within their original droplets until the latter evaporate — breakup 
mesons. Fig. |2| shows the proportion of breakup mesons for the four data sets. 
This fraction drops with increasing p c ±0 and with increasing dT/d£, because 
each of these helps the charmed quark break a flux tube and escape; the 
fraction also drops with increasing ro because a larger droplet has a longer 
lifetime and hence affords greater opportunity for escape before breakup. 
Note that the fraction of breakup mesons lies between 20% and 75%, so the 
thermalization of charmed quarks in the plasma is a very important effect. 

In Fig. § we see that breakup mesons are strongly thermalized, that is, 
their p±_ is independent of p c ±0 . (Fig. |3] shows results of a simulation with Dn 
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Figure 2: The proportion of c quarks which are trapped inside their original 
droplets until their breakup, as a function of initial transverse momentum. 
The four curves are for different values of and dV/di as in Fig. [I]. 
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Figure 3: Comparison of two populations of D mesons. The RMS p± 
of D mesons is plotted against the initial p± of the c quarks, for breakup 
mesons (dotted curve) and for fragmentation mesons (dashed curve). Here 
T* = 150 MeV, r = 1 fm, dT/d£ = 2.5 mi" 2 . 
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Figure 4: The proportion of c quarks which undergo re-ionization, as a 
function of initial transverse momentum. The four curves are for different 
values of ro and dT/d£ as in Fig. [I]. 

scattering turned off, so that hydrodynamic effects are weak and p 1 ^ is close 
to its thermal value.) The p± of the fragmentation mesons, however, does 
depend on p c ±0 and for low values of p c ±0 it in fact lies below the thermal value; 
this is due to our static model for the droplet wall, which does not include 
thermal fluctuations in its position and velocity caused by coupling to the 
plasma inside. Thus for low values of p c ±0 the proportion of fragmentation 
mesons affects p 1 ^ strongly. Having noted the dependence of this proportion 
on the parameters r and dT/di, we now understand the trends in Fig. |I[ 

These trends are strengthened by the phenomenon of re-ionization, that 
is, reabsorption by a droplet after initial hadronization. The proportion of 
quarks which undergo re-ionization is shown in Fig. |4j. This process sharply 
lowers the momentum of c quarks. The lower the quark momentum, the 
less the chances of fragmentation and the greater the probability that the 
eventual meson will emerge only at breakup of the droplet. This process thus 
converts fragmentation mesons into a second generation of breakup mesons. 
(The latter are not represented in Figs. ||] and ||.) The dependence on the 
parameters is explained as follows. Raising dT/d£ will lead to earlier release of 
the D meson, which puts it into an environment with higher droplet density, 
which will increase the chances of recapture; lowering r will increase the 
number of droplets as r^~ 3 , which again will make re-ionization more likely 
despite the decreased cross section for capture. 
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Figure 5: RMS transverse momentum of D meson vs. initial transverse 
momentum of c quark, with scattering suppressed. The four curves are 
for different values of tq and dT / di as in Fig. [I]. 

Finally, to show the importance of Dir scattering, we show results (Fig. |^) 
of simulations in which it is omitted. Comparison of Fig. |5| with Fig. [I] shows 
that Dn scattering causes a systematic increase in p^ of about 60 MeV. We 
summarize our findings at T* = 150 MeV as follows: (1) Break-up mesons 
emerge from the droplets with thermal momenta; (2) fragmentation mesons 
emerge with lower than thermal momenta; (3) Dir scattering adds momen- 
tum so that the resulting p^ is in the neighborhood of its thermal value. 

More interesting than this detailed analysis, however, is a comparison to 
the case of T* = 200 MeV, to which we now proceed. Fig. § shows that the 
RMS p 1 ^ here behaves much the same as for T* = 150 MeV, but with a sys- 
tematic shift upward of 120 MeV for given r and dT /di. (The thermal value 
of this quantity is 980 MeV, an increase of 160 MeV from T* = 150 MeV.) As 
may be seen in Fig. [7], the bulk of this effect is due to Dn scattering: Com- 
parison of Fig. [7| with Fig. || shows, in the absence of Dtt scattering, almost 
no difference between the two temperatures. Examination of the counterpart 
to Fig. [| shows (see Fig. ||) that the number of breakup mesons has decreased 
substantially, mainly because thermal motion of c quarks in the droplets is 
more effective at pushing them over fragmentation thresholds. Having seen 
that breakup mesons are thermalized while fragmentation mesons are not, 
we understand why the D mesons, upon emerging from the droplets, are no 
hotter in spite of the increase in the ambient temperature. It is only the in- 
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Figure 6: As in Fig. |T], but for T* = 200 MeV. 
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Figure 7: As in Fig. || (T* = 200 MeV), but with Dtt scattering suppressed. 
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Figure 8: As in Fig. [| but at T* = 200 MeV. 

teraction with the hot pion gas that heats the D mesons to thermal momenta 
and restores their thermometric value. 

An experimental measurement of the pj_ distribution of D mesons will 
reflect a folding of one of the curves in Fig. p] (or |^) with a distribution 
for p\ Q . As a source of uncertainty, the weak p c ±0 dependence is not very 
troublesome. More significant are the differences among the curves with 
different parameters r and dT/di. Nevertheless, as seen in Fig. 9, there is 
almost no overlap between the range of for T* = 150 MeV and that for 
T* = 200 MeV. A measurement of p^ will still provide a rough thermometer 
of the transition temperature, one that can become finer as the assumptions 
are narrowed down. 

Of the parameters we vary, we find that the initial droplet radius af- 
fects the p® most strongly; unfortunately, this parameter is the one which is 
furthest from quantitative understanding. Improvement of our results must 
therefore await a more detailed model of droplet formation in the phase tran- 
sition. 
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